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We demonstrate that under certain conditions, fractional Talbot revivals can occur in 
heterostructures of composite metamaterials, such as multilayer positive and negative index media, 
\^ ' metallodielectric stacks, and one-dimensional dielectric photonic crystals. Most importantly, 

, without using the paraxial approximation we obtain Talbot images for the feature sizes of transverse 

■ patterns smaller than the illumination wavelength. A general expression for the Talbot distance 

' in such structures is derived, and the conditions favorable for observing Talbot effects in layered 

, , heterostructures is discussed. 

^ ■ 

I—) . PACS numbers: 42.70.Qs, 42.25.Bs, 73.21.Ac 

I The Talbot effect, or the repeated self-imaging of periodic patterns, is one of the most basic phenomena in optics. 

^ ' As a consequence of Fresnel diffraction, periodic patterns can reappear (also called a revival image), upon propagation, 
at integer multiples of the so-called Talbot distance, Zt = 2a^/X, where a is the spatial period of the pattern, and A 
is the wavelength of the incident light. In addition to full revivals, fractional revivals with modified periods occur at 
' distances that are rational fractions of the Talbot distance, i.e., z/Zt = p/q, whereo and q are coprime integers. This 
Q , remarkable self-imaging phenomenon has also been demonstrated in temporal and spectralQ domains in maiw 
areas of physics from classical optics to quantum matter waves, such as gradient-index lenses^, waveguide arraysjg, 
atomic and molecular wave packets |ll Q, and Bose-Einstein condensates^. Recent investigations have revealed that 
the Talbot effect is far more than a mere optical curiosity. It is deeply connected to classical number theory and 
the intricate structure of physics Talbot effects have possible applications in optical free-space interconnects 
integrated optics, and integer factorization schemes'^ in optical computing. In parallel, several new types of composite 
CL( metamaterials with sophisticated electromagnetic properties have been developed[l3-[ll. These materials may lay 
the foundation for futuristic integrated optics and electronics. The ability to control light in certain metamaterials 
_^ has thus both scientific and practical importance. It is therefore interesting to explore Talbot effects in composite 
^ , metamaterials and potential applications involving the fractional Talbot effects, e.g., imaging in material fabrication 
fT^ ' to increase the spatial frequencies of periodic patterns. One way to control wave diffraction is to construct multilayer 
] photonic structures. Although ID periodic structures have been studied profusely, such systems still continue to reveal 
^5 ■ some interesting results jT5-|23|. Previous works have shown that photonic bands can not only control the transmission 
frequencies, but also deeply affect wave diffraction|0|-j22, leading to nondiffracting beams'2^ and high-resolution 
imaging|2l|. Typically Talbot images have been demonstrated in homogeneous media with feature sizes larger than 
' the wavelength. The higher order diffraction destroys the Talbot images when the feature sizes smaller than the 
, wavelength. In this letter, we show that through photonic engineering of wave diffraction, the fractional Talbot effect 
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^ ' can be observed for subwavelength features of transverse patterns in inhomogeneous media, such as multilayer positive 
'U^ ' negative index (PNI) materials, metallodielectric (MD) nanofilms, and low-dimensional double-dielectric (DD) 

photonic crystals. A general expression for the Talbot distance in layered heterostructures is derived and confirmed 
numerically. 

i- In any homogeneous medium, the monochromatic electromagnetic field vectors and the corresponding wavevector 

' always form a right-handed set (E,B,k) independent of the constitutive relations of the materials. Therefore, the 
• ^ continuity of the tangential field components E|| and Hy at the interface leads to the reversed sign of the normal 
component (kj^) of the wavevector when the light is incident from a right-handed medium (RHM, e > 0, > 0) to a 
H left-handed medium (LHM, e < 0, /x < 0), and vice versa, while the continuity of the tangential E|| and the normal 
- - Bj^ leads to the continuity of the tangential component (ky) of the wavevector at the interface. 

In modeling heterostructures, we assume the permittivity and permeability to be constant for the positive index 
material and frequency-dependent for the negative index material. The metal is assumed to have constant permittivity 
and frequency-dependent permeability. The frequency-dependent permittivity and permeability are given by the 



Drude model: 



w2 



^2 = 1-^-^ , t^2^l-^-^ , (1) 

where uje and are, respectively, the effective electric and magnetic plasma frequencies, and 7e and 7™ are the 
corresponding damping factors. The layer thicknesses are di and (i2, respectively, for the positive and negative index 
materials. The spatial period is d = di -1-^2, thus, e{z-\-d) = e{z) and fj,{z + d) = fJ-{z). In our layered heterostructures. 
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the index 1 refers to the positive index material while the index 2 refers to the negative index material, metal, or 
dielectric with higher refractive index. We consider TM modes {Hz = 0) and solve for the magnetic field out of 
convenience. A similar result can be carried out for TE modes. The vector wave equation for a monochromatic field 
H — H(r) exp{—iujt) is V x {e~^V x H} = (ti>/c)^/iH. Thus, the equation for the x and y magnetic field components 
becomes jl^J 

where k'j_ = k'^ + ky is the transverse wave number and k^Hx + kyHy = 0. The eigenmodes of the periodic structure 
are Bloch waves and the dispersion relation is given by 23] 

cos(/3c?) = cosh(Q!idi) cosh(a2'i2) + "^^^^ '^^^^ sinh(aidi) sinh(Q!2C?2) , (3) 

2aia2eie2 

where of — k^ — [lo / cYti^i, (i = 1, 2) and the /3 is the Bloch wave number. The existence of Bloch modes requires 
that |cos(/3c?)| < 1. It is well known that when this condition holds, the Bloch modes represent propagating waves. 
Moreover, all-evanescent modes {of > for i = 1 and 2) can also exist . In such scenarios, the Bloch modes 
represent the transmission of coupled evanescent waves. Equation ^ also represents the diffraction relation in the 
multilayer medium from which the diffraction curve, f3 versus k±, can be derived. Thus, wave diffraction in the 
multilayer structure can be constructed from a superposition of Bloch modes: 



Hx,y{v) = J dkxdkyH.j;,y{kx,ky)u*p{O)u0{z)ex.p{ikxX + ikyy + il3z) , (4) 

where Hx^y{kx,ky) is the spectrum at the z = plane and |m;3(0)P = 1. The superposition in Eq. 10} is limited 
to a single photonic band. Since the Talbot distance (Zt) is much larger than the period (d) of the multilayer 
structure, one can take the image at the distance of an integer multiple of the period closest to the Talbot distance, 
i.e. Zt = md + J « md, where m is a large integer and S <C Zt- Using the fact that up{'md) = M/3(0), Eq. Q can be 
simplified asp^l 

Hx,yir) = j dkxdkyHx,yikx,ky)exp{ikxX + ikyy + iPz) . (5) 

To produce the fractional Talbot effects from Eq. Q requires a quadratic diffraction upon propagation. This can be 
obtained by properly choosing parameters that satisfy Eq. (jSJ- 

To effectively choose the system parameters and for comparison purposes, we also evaluate the dispersion 
relation (Eq. Q) within the paraxial approximation (whereby k± <C /3o: where Pq is the center Bloch wave number 
/3o = P{k± =0)). We thus take the Taylor expansion of Eq. ^ a.t (3 = /3o- 

where Y is the right hand side of Eq. (|3J) and d is the period of the structure. In the above expansion, the odd 
derivatives vanish since Y is an even function. The fourth derivatives arising in the higher order terms are given by 



dkj 



3cos(/3od) ( d^Y\ 1 d^Y 



j=x,y, (7) 



dsin3(/3od) \ dk] J dsin(/3od) dkj 

cos(/3od) d^Yd^Y 1 d'^Y 

'dsin^iPod) dkl dkl ~ dsm{(3od) dk^dk^ ' ^' 



Near the band edges, sin(/3od) « 0, and subsequently the first terms in Eq. and Eq. JHJ grow rapidly. Thus, 
quadratic diffraction is achievable around the middle of the transmission bands with a reasonable spatial bandwidth. 
Figure n shows the diffraction curves obtained from the exact formula Eq. (|3Jl and from the quadratic approximation, 
Eq. lO, for the MD (Fig. ^) and DD (Fig. ^) stacks. In Fig. the diffraction curve in the two heterostructures 
is nearly quadratic. For comparison purposes, the diffraction curve of free space and its quadratic approximation are 
shown. All the curves are terminated at the cut-off frequency of the corresponding medium (beyond that the spatial 
components cannot be transmitted). In the PNI stack, the phase compensation effect of the negative index material 
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leads to diffraction compensation. Hence, quadratic diffraction can also appear near band edges when approaching 
the non-diffraction limit. Figure |21 shows the diffraction curves near the middle of the band (Fig. and near the 
band edge (Fig. in the PNl stack. The coincidence of the exact diffraction with the quadratic curve indicates 
quadratic diffraction inside the PNI stack. Also the diffraction in the PNl stack is much less than that of free space. 
Notice that in Figs. and 12 the cut-off spatial frequency in heterostructures is higher than that in free space. This 
property will lead to a higher resolution when transmitting images through the metamaterials. For simplicity, we 
assume the periodic pattern has the same spatial period in the x and y directions. Substituting Eq. ® into Eq. Q 
and changing the integration into double summations since the spectrum is discrete. We obtained the Talbot distance 
in the layered heterostructures as 

_ a^dsm{/3od) 

where a is the spatial period of the pattern, and y{ki, ^2) is given by 

q2y\ fQ2Y 



The coefficients read, 



A ^ [y- + sin(fcirfi + ^2^2), 



, ^1 k2 
Q = — 2 sin(fci(ii) sin(fc2d2), 



Pi 



A\ F J' fc2ei' 



4fcifc2 V fllfi2 £1^2 

Qi = ^— !- cos(A;i(ii) sin(fc2C?2) + cos(fc2d2) sin(A;ic?i), (11) 
fci k2 

where fc^ = {uj/c)'^eifii, i = 1,2. 

To confirm our theoretical predictions, numerical simulations of fractional Talbot images are demonstrated in 
Fig. O for the PNl stack, in Fig. ^ for the MD stack, and in Fig. |Slfor the DD stack. In those figures, the left plot is 
the original 2D square array while the right plot is the corresponding image at the fractional Talbot distance. The 
numerical results were obtained from Eq. jSJl where the (3 was solved numerically from Eq. not from Eq. J^. 
Unlike in free space where the paraxial approximation is required to observe Talbot images, upon carrying out the 
integration in Eq. |(SJ) no paraxial approximation was used. Equation © was used only for providing insight on how to 
effectively choose simulation parameters and for comparing the exact and paraxial results. Further, in our simulations 
the sizes of the squares in the patterns are less than the illumination wavelength. 

The interesting connection between the fractional Talbot effects and number theory can be used to explain the 
image patterns in Figs. EHHl The fractional Talbot image can be represented as a finite sum of spatially shifted 
subsidiary waves of the source field|7Ll25|: 

\ '-1 



E(x,y,^Zr) = J] (x - ^, y - o) , (12) 

^ / n 



where Z — (7/2 if q is a multiple of 4, and I — q otherwise, and 



(-1 

exp [ —1271—71^ — i2im- ) . (13) 



n=0 



For example in Fig.|3i where p/q— 1/3, the image is composed of three subsidiary waves E{x, y), E{x — a/3,y — a/3), 
and E{x — 2a/3, y — 2a/3). Thus, the spatial frequency is three times the original frequency. In Fig.|2|3 where p/q — 1/6, 
the non-zero components are &i, 63, and 65, so the image is also composed of three subsidiary waves, E{x—a/6, y— a/6), 
E{x — a/2, y — a/2), and E{x — 5a/6, y — 5a/6). Hence, the image in Fig. |3Jd has the same spatial frequency as that 
in Fig. but is spatially shifted by a half period. In realistic nanoplasmonic structures, material loss cannot be 
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avoided. The effect of material loss can be mitigated by introducing gain inside the medium'26f|. Typically, the 
damping factor of metals is much smaller than the plasma frequency J^J, 7 ^ O.Oluj^. As a demonstration, in Fig. 0| 
we compare the fractional Talbot images in the MD stack when 7 = O.Ofs"^, 7 — O.lfs^^ (loss), and 7 = — O.lfs"^ 
(gain). In the presence of loss or gain, the Talbot distance is a complex number, thus we redefine the Talbot distance 
as zt — |ztP/3?(2:t)- The Talbot image is slightly blurred when there is a loss (Fig. 0Jd), and is slightly sharper 
when there is a gain (Fig. 2t)- It is well known that at a half Talbot distance the Talbot image has a reversed 
contrast compared to the original pattern. This phenomenon is illustrated in Fig.jSjin the DD structure. For practical 
applications, we found when the error of layer thickness is less than 1%, the Talbot images can still be observed. Any 
interlayer width variation is expected to have a minimal effect since Zt ^ d. Since femtosecond lasers are widely 
used in material fabrication, for the illumination wavelength used in the DD and MD stacks, the Talbot images 
are sustainable with a wavelength variation of 7% for the DD stack and 3% for the MD stack. This approximately 
corresponds to a 100 fs laser pulses. 

In conclusion we have demonstrated subwavelength-scale fractional Talbot effects in layered heterostructures 
of metamaterials without using the paraxial approximation. A general expression of the Talbot distance in such 
structures was obtained. This expression can be used in potential applications involving Talbot effects and multilayer 
structures of metamaterials. The fractional Talbot effect can be explored in material fabrication to increase the 
spatial frequencies of the periodic patterns. This work is supported by NAVAIR's In-House Laboratory Independent 
Research (ILIR) program sponsored by the Office of Naval Research. 
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(a) Diffraction Curves 



(b) Diffraction Curves 
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FIG. 1; (Color online) Diffraction curves in the metallodielectric (a) and double-dielectric (b) stacks. The blue solid curve 
is exact from Eq. 0. The green dashed curve is the quadratic approximation Eq. The parameters in (a): A = 632 nm, 
ei — 2.66, /^i = ^2 = l,7e = 0.0 fs~^, and uje ~ 9.8 fs~^ which gives £2 ~ —9.797 from the Drude model, di = 120 nm, and 
d2 = 30 nm; in (b): A — 1550 nm, ei = 1,/ii = /i2 = l,e2 ~ 7.6, di — 60 nm and d2 = 160 nm. As a comparison, the free-space 
diffraction is shown in the red dash-dotted curve (exact) and the cyan dotted curve (quadratic). 
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(a) Diffraction Curves 
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FIG. 2: (Color online) Diffraction curves in the multilayer positive and negative index materials, ei = 2.66, /ii — l,uie ~ 
780 GHz, ujm ~ 0.8tJe, and 7e = 7m ~ 0.0 fs~^. (a) near the middle of the band, A = 5.49 mm, £2 ~ — 4.165, /12 = —2.306, di = 
0.2 mm, d2 = 0.6 mm. (b) near the band edge, A = 4.83 mm, £2 ~ — 3, /i2 = — 1.56, di — d2 = 0.2 mm. The coincidence of the 
blue solid curve (exact) with the green dashed curve (quadratic) indicates the diffraction is quadratic in the PNI stack. As a 
comparison, the free-space diffraction is shown in the red dash-dotted curve (exact) and the cyan dotted curve (quadratic). 
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Object: 2D Arrays 



(a) Fractional (1/3) Talbot Image 
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(b) Fractional (1/BJ Talbot Image 
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FIG. 3: (Color online) Fractional Talbot images in the PNI stack. The size of the square is 2 mm. (a) z — Zt/S, Zt = 1.2 m, 
all parameters are the same as those in Fig. |2K (the middle of the band), (b) z = Zt/G, Zt ~ 7.3m, all parameters are the 
same as those in Fig. (the edge of the band). Due to the phase compensation of the negative index medium, the Talbot 
distance is long. 



Object: 2D Arrays 
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FIG. 4: (Color online) Fractional Talbot image in the MD stack in the presence of loss or gain. The size of the square is 0.4 /im, 
z^= Zt/8, and = 9.8fs-\ (a) 7 = O.Ofs"\ Zt = 125.1 Atm; (b) 7 = 0.1 fs"^ (loss), Zt = 125.3 Mm; (c) 7 = -0.1 fs"^ (gain), 
Zt = 125.3 /im. All other parameters are the same as those in Fig.^. 
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Object 2D Array 
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FIG. 5: (Color online) Fractional Talbot image in the DD stack. At a half Talbot distance the image has undergone a contrast 
reversal compared to the original pattern. The size of the square is 1/im. z — Zt/^, Zt ~ 12.9 /im, all other parameters are 
the same as those in Fig.^. 



